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Presence of a (1x1) oxygen overlayer on bare ZnO(OOOl) surfaces and at Schottky 
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The atomic surface and interface structure of bare and metal-coated ZnO(OOOl) Zn-polar wafers 
were investigated via surface x-ray diffraction. All bare samples showed the presence of a (1x1) 
overlayer of oxygen atoms located at the on-top position above the terminating Zn atom, a structure 
predicted to be unstable by several density functional theory calculations. The same oxygen overlayer 
is clearly seen at the interface of ZnO with both elemental and oxidized metal contact layers. No 
significant atomic relaxations are observed at surfaces and interfaces processed under typical device 
fabrication conditions. 
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ZnO has many important technological applications 
in fields as diverse as catalysis, gas sensing, corrosion 
prevention, and optoelectronics.^ A precise knowledge 
of the structural properties of the surfaces and inter- 
faces involved in these applications is a prerequisite for 
their control and optimization. However, despite a large 
number of experimental and theoretical studies, the sur- 
face and interface structure of ZnO remain a topic of 
intense debate, with many alternative structures having 
been describedi~— This lack of consensus arises primar- 
ily from the competition of various surface stabilization 
mechanisms and a host of different preparation proce- 
dures and ambient conditions, resulting in a complex 
phase diagram. ^^ 

In this letter, we report on a systematic study of the 
atomic structure of bare ZnO surfaces and metal Schot- 
tky contacts to ZnO prepared under typical device fab- 
rication conditions. Schottky contacts are important 
building blocks in many electronic devices, and an under- 
standing of their interface structure is important, since 
electronic and structural properties are usually strongly 
correlated. This is particularly true for semiconductors 
characterized by highly ionic bonding (e.g., ZnO), in 
which small atomic displacements can potentially result 
in large changes in electronic behavior. 

Single-crystal ZnO (0001) wafers were hydrothermally 
grown along the +c-axis by Tokyo Denpa Co. Ltd. and 
epi-polished to low miscut angles (< 0.1°).^^ For mea- 
surements on bare surfaces, samples no. 1-3 were ul- 
trasonically cleaned using organic solvents, dried in N2, 
and subsequently exposed to air. Polycrystalline Schot- 
tky contacts were deposited on other samples after the 
same cleaning procedure, either as plain gold by ther- 
mal evaporation (samples no. 4,5) or non-stoichiometric 
iridium oxide (IrO^;) layers by eclipse pulsed laser de- 
position in an oxygen ambient (samples no. 6,7). The 
latter method has been shown to consistently produce 
high quality Schottky contacts to ZnO.^^ 

Surface x-ray diffraction (SXRD) offers the unique ca- 
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FIG. 1. (color online) Atomic structure of a bulk- like 
ZnO (0001) surface, (a) top view with the surface unit cell. 
(b) side view along the [110] direction. Note the presence of 
two ZnO double layers (DLs) within the unit-cell height c. 



pability to access the buried interface structures of the 
Schottky contacts non-destructively, and it provides pi- 
cometer accuracy for the determination of atomic po- 
sitions. SXRD experiments were carried out on the 
Materials Science beamline X04SA at the Swiss Light 
Source (SLS) and on the GeoSoilEnviroCARS beamline, 
13-BMC, at the Advanced Photon Source (APS) using 
12.398 keV and 15.000 keV photons, respectively. A 
PILATUS lOOK pixel detector^^as ^^^ ^^^^ f^^ f^g^ ^^^ 

reliable data acquisition. 

In-plane line scans along high-symmetry directions 
gave no evidence for surface or interface reconstructions 
on any sample. The out-of-plane crystal truncation rod 
(CTR) measurements revealed a 6-fold rotational sym- 
metry of the diffraction pattern {p6mm)^ even though a 
perfectly bulk-truncated ZnO crystal is expected to have 
p3ml symmetry. However, there are two ZnO double 
layers (DLs) in the ZnO unit cell (Fig. [T]), both of which 
result in a chemically identical termination layer but a 
180° rotation of the p3ml diffraction pattern. The overall 
p6mm symmetry is therefore consistent with the presence 
of 1/2 unit-cell height terrace steps and equal propor- 



tions of both terminations on the surface. Atomic force 
microscopy measurements on hydrothermal ZnO(OOOl) 
Zn-polar surfaces cleaned in the same way have previ- 
ously showed evidence of triangular islands and pits with 
180° rotation between triangles on terraces separated by 
a single DL step.— 

For each sample, at least 8 symmetrically inequiva- 
lent {p6mm) CTRs plus several equivalent rods were 
recorded, typically resulting in 400-700 averaged struc- 
ture factors (SFs) per data set and systematic errors 
of 5-10% between symmetry equivalents. A specialized 
module to model SXRD data for the genetic algorithm 
refinement program GenX^^ was used to fit the mea- 
sured SFs. Utilizing differential evolution algorithms, 
GenX efficiently avoids becoming trapped in local min- 
ima. To verify the reproducibility and uniqueness of a 
solution, all fits were repeated at least ten times with ran- 
dom parameter initializations. To give the low-intensity 
regions of the CTRs a similar weight in the fit as the 
high-intensity points close to the Bragg peaks, a loga- 
rithmic i?- factor (i^iog) was employed for the fitting fig- 
ure of merit (FOM).^^ All final fit results are also given 
in terms of the standard crystallographic i?- factor j^ 

To identify the characteristic structural features, a 
large number of different models with varying degrees 
of complexity based on the bulk-like surface struc- 
ture shown in Fig. [1] were tested initially. Atomic 
z-displacements (Az), occupations (pocc), and Debye- 
Waller (DW) factors of all atoms in up to four ZnO DLs 
were allowed to vary freely. The local pSml surface sym- 
metry only permits atomic movements Az along the z- 
direction and allows atoms to be present in just three 
different positions in the unit cell: at (x, y) = (0, 0) (fee 
hollow), (1/3,2/3) (on-top), and (2/3,1/3) (hep hollow). 
In addition to a completely bulk-like structure, we there- 
fore tested various models with adatoms (both oxygen 
or metal atoms) located at some or all of these allowed 
locations and free to move within one unit cell away from 
the surface^. 

The results of these survey fits draw a consistent pic- 
ture of the surface or interface structure across all sam- 
ples, including those with metal layers. The occupa- 
tion parameters of Zn reveal a sharp surface/interface 
with only one partially occupied atomic layer, where 
j9occ(Zn) < 40%. The corresponding oxygen occupation 
within this DL is, however, close to unity, as is the Zn 
occupation one DL below. Additionally, we usually ob- 
serve a comparable partial occupation of oxygen atoms 
in the next DL above the incomplete Zn-layer in the on- 
top position, while the occupations of the corresponding 
Zn atom and any adatoms in fee and hep hollow sites 
are negligible in that layer. Hence, there is strong ev- 
idence for an oxygen overlayer on-top of the nominally 
Zn-terminated surface, while there are no indications for 
the presence of any other ordered adsorbed species or, in 
particular, the ordering of metal atoms in those samples 
with metal Schottky contacts. Apart from the oxygen 
overlayer atoms, which show an outward relaxation of ap- 




FIG. 2. (color online) Measured diffraction data (open sym- 
bols) and corresponding calculated intensities (lines) for two 
representative CTRs on all samples. Error bars are smaller 
than the data points and have been omitted. 



TABLE L Final R\og and R values on all samples [xlO ^]. 



Sample No. 
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3.87 
6.22 


1.54 
6.37 


1.42 
6.71 


1.75 
7.56 


1.41 
4.67 


2.03 

8.25 


1.49 
4.11 



proximately 50 pm, there are no significant movements 
(Az < 10pm), even in the topmost atomic layers^*^*^ The 
DW factors converge to values close to those reported for 
bulk crystals ^^ and manually modifying them has no sig- 
nificant infiuence on the optimized values of any other 
parameters. 

For a direct comparison of all samples, the trends iden- 
tified in the survey fits were captured in one model, which 
was systematically applied to all data sets. It consists of 
a bulk-like ZnO unit cell with an extra oxygen overlayer 
added on top of the terminating Zn layer, and includes 
occupation parameters, DW factors, and z-displacements 
for the two topmost DLs (Znl, 01, Zn2, 02; see Fig. [H 
and the overlayer (Oql)- This model yields good fits on 
both bare and metal coated surfaces. Fig. [2] shows the 
excellent agreement between measured diffraction data 
and simulated CTR profiles for all samples, and the cor- 
responding final FOM values are summarized in Table HI 

Fitted z-positions and occupation parameters are plot- 
ted in Fig. [3l The presence of the oxygen overlayer is 
clearly seen, as the occupation of the two topmost oxygen 
layers (Oql, 01) is equal to those of the corresponding 
zinc layers immediately below (Znl, Zn2). The low occu- 
pation (< 35%) of Znl and Oql is most likely attributed 
to the presence of small islands on an otherwise atomi- 
cally fiat surface. Thus, a large fraction (> 65%) of the 
01 atoms form the overlayer on top of the fully occupied 
Zn2 layer where the islands are absent. No significant 
atomic movements within the ZnO are observed (Zn2, 
02: Az < 5pm), but the tendency for outward displace- 
ments of the 01 and Oql atoms suggests an overlayer 
bond length larger than the Zn-0 bond length in the 
bulk.i^ 

The (1x1) oxygen overlayer on the bare Zn-polar ZnO 
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FIG. 3. (color online) Fitted z-positions and site occupations 
for all samples. Average error bars for the parameter values 
are shown next to the atom labels. These are defined by a 
5% increase in the FOM value. Dotted vertical lines mark the 
nominal bulk positions. 



surface is most likely in the form of hydroxyl groups 
(OH), previously observed by x-ray photoemission spec- 
troscopy (XPS),^^^^^ although we lack sufficient sensitiv- 
ity to detect hydrogen with SXRD. XPS measurements 
on the same hydrothermal ZnO material indicate a sta- 
ble OH coverage of at least 1 monolayer (ML) on the 
Zn-polar face^ However, the formation of this fully oc- 
cupied (1x1) overlayer is not well understood, as such 
an arrangement does not obey the electron counting rule, 
and is predicted to be thermodynamically unfavorable by 
density- functional theory (DFT).^^^^^ Furthermore, this 
fully occupied overlayer does not appear in any of the 
calculated phase diagrams for ZnO, which predict no sta- 
ble phases with greater than 1/2 ML OH coverage and 
OH adsorption at fee- hollow rather than on-top sites i~— 
Our results therefore indicate that alternative stabiliza- 
tion mechanisms may play a role at surfaces prepared 
under typical device fabrication conditions. 

Interestingly, the ordered (1x1) oxygen layer remains 
intact when the surfaces are covered with plain or oxi- 



dized metal Schottky contacts. This is consistent with 
ab initio calculations of the adsorption of Cu atoms on 
polar ZnO surfaces and the role of chemical bonding 
at metal- ZnO interfaces, which indicate that metal-zinc 
bonding is unfavorable and associated with ohmic rather 
than Schottky contact behavior. ^^^^^ The question as to 
whether metal adsorption occurs on top of the surface 
OH groups or by replacing the H atoms remains open, 
as the detection of hydrogen at buried interfaces is ex- 
tremely challenging. 

In conclusion, SXRD showed that the atomic structure 
of the bare ZnO (0001) Zn-polar surface prepared under 
typical device fabrication conditions has a bulk-like ter- 
mination with no significant atomic relaxations. Most 
interestingly, a (1x1) overlayer of oxygen atoms on top 
of the terminating zinc atoms was observed, consistent 
with XPS measurements, but at odds with DFT calcu- 
lations. At bare surfaces, this (1x1) oxygen overlayer 
is most likely associated with the presence of hydroxyl 
(OH) groups. Significantly, the (1x1) oxygen overlayer 
remains in place following the fabrication of both plain 
and oxidized metal Schottky contacts. 
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